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In Part I of this paper we formulated the problem of error detection with quantum 
codes on the completely depolarized channel and gave an expression for the probability 
of undetected error via the weight enumerators of the code. In this part we show that 
there exist quantum codes whose probability of undetected error falls exponentially with 
the length of the code and derive bounds on this exponent. The lower (existence) bound is 
proved for stabilizer codes by the counting argument for classical self-orthogonal quaternary 
codes. Upper bounds are proved by linear programming. First we formulate two linear 
programming problems that are convenient for the analysis of specific short codes. Next we 
give a relaxed formulation of the problem in terms of optimization on the cone of polynomials 
in the Krawtchouk basis. We present two general solutions of the problem. Together they 
give an upper bound on the exponent of undetected error. The upper and lower asymptotic 
bounds coincide for a certain interval of code rates close to 1. 
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Index Terms — Probability of undetected error, self-orthogonal codes, polynomial method. 



1 Introduction 

In part I of this paper we defined the undetected error event for transmission with quantum 
codes over a completely depolarized channel and explained a way to compute its probability 
via the weight enumerators of the code. This part is independent of part I once we agree on the 
definitions of a quantum code, the channel and error correction, and the weight enumerators. 
The main results of part I can be summarized as follows. Let Q be an ((n, K)) quantum code, 
i.e., a linear K-dimensional subspace of TC n := (C 2 )®". Let 

n 

B(x,y) = X>^ n -y 

i=0 



B ± (x,y) = £>, 



V- 



=0 



*Los Alamos National Laboratory, Group CIC-3, Mail Stop P990, Los Alamos, NM 87545. 
^Bell Laboratories, Lucent Technologies, 600 Mountain Avenue 2C-375, Murray Hill, NJ 07974. 
*Los Alamos National Laboratory Group CIC-3, Mail Stop P990, Los Alamos, NM 87545. 
^Department of Electrical Engineering-Systems, Tel Aviv University, Tel Aviv 69978, Israel. 



be the (Shor-Laflamme) weight polynomials of Q, where the weight distributions Bi, < % < n, 
and Br-,0 < i < n, are given by (1.9) and (1. 10), respectively^. Then the probability of 
undetected error equals 

n 

Pue(Q, P ) = X> 4 X " Bi)(J)\l ~P)\ (1) 

i=l 

where p is the probability for a factor Tj in the error operator E = t\ ® • • • <S> r n to be nonidentical 
(this probability does not depend on i). In some situations this formula should contain an 
additional constant factor (see Theorem 1.6) which throughout this part will be omitted. 

The goal of this part of the paper is to derive bounds on P ue (Q,p) for the best possible code 
Q with given parameters. More specifically, in Part I we defined the quantity 

P ue (n,K,p) = min P ue (Q,p). 

dim(Q) = K 

We will derive upper and lower bounds on P ue (n,K,p). Just as in the classical case, this 
probability falls exponentially with n; therefore, let us also introduce the exponent of undetected 
error 

E(n,R Q ,p) = --log 2 P ue (\2 R Q n ],n,p), 
n 

where Rq = — ^ — is the code rate. Speaking of asymptotics, we are interested in upper bounds 

E(Rq,p) = limsup£'(n, Rq,p) 

ri— >oo 

and lower bounds 

E(R Q ,p) = liminf E(n,R Q ,p) 

n—>oo 

(this corresponds to lower and upper bounds on the probability P ue (n,K,p), respectively). Let 
E(Rq,p) be the common limit of these functions, provided that it exists. 

Throughout the paper 

T q (x,y) = xlog q (q- 1) - xlog q y - (1 - x) log g (l - y) 
H q (x) = T q (x,x). 

In the classical case, error detection has been extensively studied. The probability of unde- 
tected error in the classical case is defined in (1.1); its exponent E^ d '(R,p) is defined similarly 
to the above. Best known lower bounds on E^ >{R,p) (upper bounds on the probability) were 
derived in J9| building upon the Varshamov-Gilbert-type existence arguments. We consider the 
binary case only. Let R vg (x) := 1 — H2{x) be the Varshamov-Gilbert function and 5 vg {x) its 
inverse. Also let R(x) be the function giving the best known upper bound on codes |12| and 
5(x) its inverse. It is easy to prove that there exist binary linear codes with Ai < n(")2 , 

References (1.9), Theorem 1.6, and so on point to equations, theorems, etc., in the first part of this paper. 



where Ai is the number of vectors of weight i in the code. Substituting this in (1.1), we obtain 
the lower bound of ||: 

E (d) (R v) > iT 2 (5 vg (R),p) 0<R<1- H(p) 
} ~\l-R l-H 2 (p)<R<l. 

Upper bounds require more involved arguments 0, [[ll]]. The results have the form 

E ( d)(R s < f 1 - R ~ H 2 {6{R),P) + T2(6(R),p) 0<R< R{p) 
1 ' ~ \\-R R{p)<R<l- 

In this paper we derive analogous results for the quantum case. In Section g we prove the 
existence of quantum stabilizer codes with bounded above weight enumerators; a substitution 
in (|l]) yields lower bounds on E(R,p). In this part we rely upon the results of [L|] on quaternary 
self-orthogonal codes. Then we move on to lower bounds on P ue (n,K,p). In this part we 
employ the linear programming technique. In Section ^ we formulate a linear programming 
(LP) problem with the objective function P ue (n, K,p). Though in examples this problem gives 
good lower bounds (which can be found by solving it on a computer), analytically it is difficult 
to deal with. Therefore, in the second part of this section we propose a relaxation of the problem 
which enables us to derive general bounds. This part of the paper is based on an application 
of the LP approach in the quantum case [|| in conjunction with the methods of j|], [11]. The 



results include two upper bounds on E(R,p). These bounds show that for the rate Rq in a 
certain neighborhood of 1, dependent on p, the exponent E(R,p) is known exactly. For lower 
rates, the bounds are in general location, i.e., there exists a segment in which each of them 
is better than the other. This part of the paper is technically the most involved. We chose 
to formulate the results for arbitrary q (instead of concentrating on q = 4), the reasons being 
that once we look at q > 2, it does not make much of a difference whether it is 4 or anything 
else, and that this is helpful in studying error detection of nonbinary classical codes on which 
we plan to report elsewhere. Moreover, the theory of quantum codes generalizes to larger q 



[g], 1 14], though the presentation is somewhat less systematic and the results more scattered 
than for binary quantum codes. 

Some further remarks on notation. Throughout the paper F = F4 = {0, l,u;,u; 2 }. The 
Krawtchouk polynomial is given by 



1=0 W 



»:*)(«- ir 



(the implicit parameter n - the length of the code - is usually clear from the context). Properties 
of Kk{q;x) used throughout the paper are summarized in the appendix. As remarked above, 
by Rq we denote the rate of the quantum code. We also use two associated numbers |C| = 
Y^i=o-Bi an d |C I = Y^i=o^t'^ i n t ne case of stabilizer codes they are equal to the size of 
the two underlying classical codes, C and C 1 - (see Part I). Likewise, let R = (1/n) log 4 |C|, 
R = (1/n) log 4 |C |. The rate Rq and these 2 quantities are connected by the following 
relations: 

Rq = 2R 1 - -1 = 1-2R, R + R ± = 1, < R < - < R 1 - < 1. (2) 



2 Upper (existence) bounds on P ue (n,K,p) 

In this section we show that there exist quantum codes for which probability of undetected error 
falls exponentially for all rates < Rq < 1, and bound this exponent below More specifically, 
we prove the following theorem. 

Theorem 1 

E(R p)>/ r4 ( iT 4 _1 (K 1 -^))^) <R Q < 2(1 -H 4 (p))-1; 
rU} ' P '~\\{1- Rq) 2(1 - H 4 (p)) -1<R Q <1. 

To prove this theorem, we restrict our attention to quantum stabilizer codes. In analogy 
with the classical case, we show that there exist sequences of codes Q of growing length n and 
size K = 2 Q with weight distribution 

Bt<n 2 ( n )3 i 2 k Q- n . 

(in fact, n 2 can be easily replaced by n). 

Since the weight distributions B L , B^ correspond to classical quaternary code, we prove this 
estimate by considering families of quaternary self-orthogonal codes. Let C C F n . Throughout 
the section we denote by C 1 - a linear code dual to C with respect to the standard dot product 

(-,•)• Let 

S n ,k = {C C F | C even linear code, dimC = k} 
S^ k = {CdF\C L £ S n , k }, 

where k < n/2 by (||). 



Below we use the following three results from [13]. 



Lemma 2 Let C C F n be an even linear code. Then C is self- orthogonal with respect to the 
inner product a * b = Y^i=i( a i^i + a ^«)- 



Lemma 3 Let C C F n be an even linear code and C be its dual. Then the number of even- 
weight code vectors in C 1 - equals 7,(4" + (— 2) n ). 



Lemma 4 Let C C F n be an even linear code and C be its dual. I/a£ C has even (odd) 
weight then the coset C + a is formed by vectors of even (odd) weight. 

Existence of codes with bounded distance distribution will follow from the following lemma, 
based on Lemmas ||-|]. 



Lemma 5 Let v G F n be any even-weight vector. The number of codes from S containing v 
does not depend on v. 

Proof. Let us count the number of linear [n, k] codes containing v. Let C\ = (v) be the [n, 1] 
code and w an even-weight vector distinct from v such that (w, v) = 0. By Lemma all the 
cosets C\ -faw, a £ F, are even; so by adjoining w we obtain an even [n, 2] code C^- By Lemma 
this can be done in 

\{f- l + (-2)")-l 



1 
4 
ways independent of v. Similarly C2 can be extended to an even [n, 3] code in 



1 
16 



' ' ?~ 2 + (-2) n ) - 1 



ways. Continuing in this manner, we obtain all even [re, k] codes that contain v. It is obvious 
that their number does not depend on a particular choice of v. □ 

Theorem 6 The family S Ut k of even [n, k] codes contains a code C with weight distribution 
Bi(C) < n 2 Bi, 1 < i < re, where 



l^ 4 n-l + (_ 2 )nj _l \l 

is the average weight distribution of codes in S n ^- 

Proof. Let N = |<S nj fc|. By Lemma g every even vector v is contained in one and the same 
number, say L, of codes from S n ^k- So computing the total number of all vectors in codes from 
S n: k in two ways, we get 



N 



+ (IjV" 1 + (-2)") -1\L = A k N, 



or 

i(4"' 1 + (-2) nN ) -1 

N = L-± E '- . 

4 fc -l 

Let Bi(Cj) be the number of code vectors of weight i in the j-th code from S n: k- Then we have 

JV 



E*M = (")«■ 



3=1 

Hence the average over S n ± number of codewords of weight i is -4 (™)3 l = Bi, as claimed. The 
number of codes C £ 5 n) fc such that Bi(C) > n 2 Sj for a given i is not greater than 



SiW_ii( 4 +( - 2)B J- 1 _JL^ 



?? 



2 5. n 2 4 fc - 1 re 2 ' 



Hence the number of codes C such that Bi(C) < n 2 Bi for all 1 < i < n is at least 

W -£_w(i-I- 

Now let us use this result to prove that the family S^ k also contains codes whose weight 
distribution is bounded above by a polynomial factor times the average weight distribution 
in S^ k . This will enable us to prove Theorem ffl. In this part we rely on the Mac Williams 
identities. We will need the following lemma. 

Lemma 7 Let n be an even integer. Then 

n/2 



J2(fy 2l K t (A-,2 l )=2^(f)(-3y. 



Proof. By (39), the sum J2"=o (^)3 2i ^r(4; 2i)K s (4; 2i) is the coefficient of y r z s in 
n/2 . . 

S ( r ) 32i(1 + 3y ) n_2i ( 1 " y) 2 ^ 1 + 3z ) n ~ 2i ( 1 - z f l 

i=0 ^ *' 

= 5([(l + 3y)(l + 3*) + 3(l-|/)(l-*)] B 

+ [(l + 3y)(l + 3z)-3(l-y)(l-z)r) 

= i([4 + 12yzr + 2"[-l + 3(y + z)+3yzr). (3) 

It is clear that the first term in (3) contributes only to coefficients of y r z r . Consider the second 
term: 

2 n [-l + 3(y + z) + 3yz] n = 2 n ^ ( n \ 3 n ~< ' y n ~ % 'z^ 1 ; (3(y + 2) - 1)* (4) 

i=0 ^ ' 

Since we are interested only in the coefficient y°z t , in the sum (4) we put i = n. This gives 

2"(3( y +z)-ir = 2" £ ( n Vsr E V v-J 
= 2"E/'E(")Q(-3)v-' 

The coefficient of y z f in this sum equals 2 ra (™)(— 3)*. □ 

Theorem 8 In S^ k there exists a code C with Bj- < n 2 B^~, 1 <i < n, where 

B^~- ( n ) 3* (l + ((- 1 ) t2 "' 1 - 1 )( 4fc - 1 ) ^ (5) 

B% ~ ^{i) d \ 1+ i( 4 «-i + (-2)»)-l ) {) 

is the average weight distribution of codes in S^ k . 



Proof. As in Theorem [6|, it suffices to compute the average. Using the Mac Williams identities 
and the fact that \S^ k \ = |5 n ,fc|) we have, for a given t, 

N N n 

E B ^ c f) = ziEE^ow 4 ;*') 

j=l j=l j=0 

n N 



4fc 



i=0 j=l 



n/2 



i=0 



,1,1 i. . . 
= ^(iV^(4;0)+L^( 2 n 3 2 ^(4;2,)-LK t (4;0)). 

From Lemma ^ it follows that for t > 



i=i 



4*' 



|:y+I2"- J (-:]) , (") -L( n ):V 



4 fc -l 



1 + ((-l)^™- 1 - 1)- 

VV ' ^( 4 +(- 2 ) n ) 



Hence 



A-^XX^ 



i=i 



□ 



The proof is completed as in Theorem |FJ. 

Now we are in a position to prove Theorem FT]. Let Q be a quantum code satisfying Theorems 
H and ||. Note that the second term in the expression (||) for Bj- vanishes as n grows; so starting 
with some value of n the weight distribution B± is bounded above as 



£f <2n 2 . 3*2 



/? 



linko-n 



(6) 



Let us compute P U e(Q,p) for large n relying on this inequality. We have 



p\ % 



Pue(Q, P )=Y,( B t-Bi) (5 )(1-P) 



i=0 



sE^dVa-pr 



i=l 
n 

,2\ .y.-,'- 



i=l 



;l„/^2 fc «-™ +1 (^)3'f^) (l-p) r 



(7) 



The exponent of the summation term equals 



n 



Tui)-^p-(i-i)**u(i-P) + \<i-XQ) 



(8) 



where we have omitted the o(l) terms. The expression in brackets attains its maximum for 
(i/n) = p. Note also that the right-hand side in (|6|) behaves exponentially in n; the exponent 
approaches when the quotient (i/n) — > 5 vg (R- L ). Hence as long as (i/n) > 5 vg (R ± ), the sum 
in (7) asymptotically is dominated by the term with i = \np\. Thus as long as 1 — H^(p) < 
R 1 - = ^(1 + Rq), we have 

E(R Q ,p)< l -(l-R Q ), 

i.e., the second part of the theorem. Otherwise, the maximum moves outside the summation 
range, so the largest term asymptotically is the first one in the sum, i.e., the one corresponding 
to i = [n5 vg (R- L )\. This gives the first part of the theorem. □ 

Note that we have proved a stronger fact about quantum codes than the one actually in 
Theorem Q, namely, that there exist stabilizer codes both of whose weight distributions Bi and 
B± are bounded above by the "binomial" term n 2 ( n )3 l 2 k Q~ n . 



3 A linear program for quantum undetected error 

The approach leading to best known lower bounds on the probability of undetected error in the 
classical case has been the linear programming one H , [|ll|] . In this section we develop a similar 
technique for the quantum case. First we formulate two theorems that enable one to obtain 
good lower estimates on P ue (n,K,p) for finite n. Then we formulate a relaxed LP problem 
which is not as good for finite n but lends itself to asymptotic analysis. 

For the reasons outlined in the introduction we will study a general alphabet of size r. An 
r-ary quantum code Q is a ET-dimensional linear subspace of C r . Without going into details we 
say that one can associate with Q two weight distributions, Bi and B^~, < i < n, connected 
by the q-aiy Mac Williams identities, q = r 2 . Furthermore, 



1 - 



j=0 



As above, we use the notation 



it 



i=0 i+0 

R = (l/n) log, \C\, R ± = (l/n) log, \C^\; 

these numbers and the rate Rq are again related through @. We have for the probability of 
undetected error 

n 
Pue(Q, P ) = Y,( B f ~ ^(^l)^ 1 ~P) n ~ K W 

i=i q 

Our first result is given by the following theorem. 



Theorem 9 Let Q be an {{n,r nR Q)) quantum code. Let q = r 2 , R 1 - = (1 + Rq)/2. Let 
Z( x ) = J^i=o z i^i(Q'i x ) an d Y(x) = X^r=o^-^*( ( ?' 2 ') ^ e polynomials such that 



Z(j)-Y(j)+y + yj q 



nR 1 



<a nR^(JP_y ( ,_^n- J fq-i-qPV 



q-l 



(i- P y 



q-l 



, 1 < 3 < n, 



< 



Zo, yo ^0; zj > 0, yj >0,l<j<n. 



Then 



P ue (n, r nR Q,p) > q ~ nR± (z q nR± - Z(0) + Y(0) - y ) + (1 - p) n - q 



nR 1 



Proof. Using the MacWilliams identities we can rewrite (||) as follows: 

1 



Pue(Q,P) 



Y.ZJ 



ir'J-i ^ 3 



\C^\v[^\\l-vT-i 



q-l 



E^W^rJt 1 -^ 



\n—t 



t=l 



/? 



ICJ-i — \j 






p>{l-p) n - j . 



The middle term here is calculated using the generating function (39): 

n . ^ 

p \*m xn-t (q-^-qpv 



£*<«)(-^) a -*>)-' 



i=0 



q-l 



Hence 



Pue(Q,P) 






\c^\(^-Y(i- P r-j 



1 — qp\i 



1 



+ (1-P) n 



hU^ {1 ~ pr 



|CJ-| 

1 



+(l-p) n - 



ic-Hf-^Ya-pr-' 



1 — gp\ i 



1 



\CM' 



(10) 



(11) 



Thus, P U e{Q:P) is a linear form of the coefficients B^ which we have to minimize. We can 
formulate the following LP problem: 



3=1 



min < > i? 



j 



ic ^Myv^ 



1 ~~ qp\J 



q-l 



(12) 



subject to the restrictions 

Bf > 0, 1 < j < n 

n 

E B f = \ c± \ - 1 



< % < n 



J2BfK t (q;j) > -(^j( q -l)\ 1 

BJ-lC^-f^Bl-K^i) > ( n )(q-iy, l<j<n. 

The last inequality follows from Theorem I.2(i). 

Now the theorem follows by the LP duality. Indeed, the dual program has 2n + 1 variables 
(zo, z\, . . . , Zn) and (yi,y2, • • • , y n ) of which zq can take on any value and all the other variables 
are nonnegative. The dual objective function has the form 



max 



{z (i^i - 1) - f> - ^V - ir} 



subject to restrictions 



ZQ 



+ Y J {zi-y i )K i {q;j)+y j \C A 



i=\ 



Any feasible solution of the this problem gives a lower estimate of P ue (n, K,p). Let us intro- 
duce the polynomials Y(x) = ^2 r l = QyiKi(q;x) and Z(x) = ^«=o z iKi(Q'i x )- Since -Kofe^) = 1, 
this implies our claim. □ 

Sometimes it is convenient to rewrite the linear problem via the enumerator Bj instead of 
Bj~. The proof of the following theorem is similar to the above. 

Theorem 10 Suppose Q,n,r,RQ,q, Z(x),Y(x) have the same meaning as in Theorem^. Let 
R = (1 — Rq)/2. Suppose that 



Z(j)+Y(j) -y - Vi 



nR 



s( i rn JS )'-«" , '(irT)'< 1 -j>r', i<j<», (is) 



<, 



Then 



zo, Vo ^0, Zj > 0, yj > 0, 1 < j < n. 

Pue(R q ,p) > q" nR (z q nR - Z(0) - Y(0) + y ) - (1 - p) n + q'™. 

10 



Though the LP problems of Theorems |9| and 10 enable one to find bounds for short codes 
with the use of computer, they are not easy to analyze in general. The reason for this is that 
the sign of the quantities on the right-hand side of (10) or (13) alternates. This significantly 
complicates checking feasibility of a putative solution. For this reason below we take on a 
different approach which, though it does not yield optimal solutions for the LP problem, gives 
rise to good asymptotic upper bounds on E(R,p). 

Theorem 11 Let Q be an r-ary quantum code with weight enumerators Bi and B^~. Let 
h(i), 1 < i < re, be a real-valued function and 



Z{x) = S ^z i K i (q]x), (q 



r 2 ) 



i=0 
be a polynomial that for 1 < i < re satisfies the conditions 



Then 



(i) Z(i)<h(i); (14) 

(ii) Zi\C L \ - Z(i) > 0. (15) 



J2(B±-B i )h(i)>z \C ± \-Z(0). (16) 

j=0 



PROOF. The proof will follow from the following chain of relations: 

n n 

Y^Bj- - B^hd) > J2^ - B % )Z(i) 



i=0 i=0 

i=0 ' ' j=0 

n 1 n n 

= Y. B iTc\^ Kl{q ' J)Z{t) ~^ BiZ{i) 

j=0 I' i=0 i=0 

n n n n 

= E b itc\ z > -T. B ^) = J2 B i(\ c± \ z i - Z (J)) 

j=o I I i=0 j=0 

> \C x \z - Z(p), 

where the first inequality follows by (i) and the obvious Bq = Bq = 1; step (a) is implied by 
the Mac Williams identities, in (b) we use (41), and the final inequality follows by (ii) and the 
fact that Bj > 0,B = 1. □ 

If (15) is replaced by the condition 

Z(i)-\C\ Zi >0,i = l,...,n, (17) 



11 



then by a similar argument one can prove that 

n 

Y,{Bi - B l )h(i) > Z(0) - \C\z . 

i=0 

We wish to stress the difference between the conditions on Z(x) in this theorem and in the 
classical (non-quantum) case ||. The standard condition Zj > 0, 1 < j < n, is not needed 
to prove (|T6|); it is replaced by related though not equivalent conditions (15), fli~7|). In the 
situation when one bounds above the size of the quantum code, the corresponding inequality is 
\C\ < Taax.i<j< n (Z(J)/zj); see ||] for details. 



In the next section we use Theorem 11 to derive asymptotic bounds on E(R,p). 



4 Lower Bounds on P ue (n, K,p) 

In this section we prove two lower bounds on the probability of undetected error that are valid 
for any quantum code of given length and size. The bounds are derived by a suitable choice 
of polynomials in Theorem [TJ]. The results are similar in spirit to ||, P], 111]]. In this section 
J = q- 1, 

h(x) = (p/ 1 ni-p) n - x , (18) 

and as usual < p < y/q. 

4.1 An Aaltonen-MRRW-type bound 

In this part we rely on the technique in JI^q extended to arbitrary q in [U, and apply it in a 
way similar to ||, ||. Let 

t (z):=?--^±z--^/ 7 z(1-z) (0 < z < (j/q)). (19) 

By |1^| , m Rq(5) := H q (To(5)) is the maximal asymptotically attainable rate of a g-ary (clas- 
sical) code with relative distance 5. This is proved by studying Delsarte's linear programming 
problem with the polynomial 

f t {x) = —^—AK t {q;x)+K t+1 {q-x))\ 
(a — x) 

where t = [tit\, t = r (5), and a is the smallest zero of K t (q;x) + K t+ i(q;x). Conversely, the 
function 5q(R) := tq(HZ' 1 (R)) gives an asymptotic upper bound on the minimum distance of 
g-ary codes of rate R. (An interesting remark is that the function to(-) is involutive; this is 
intimately related to the self-duality of the Hamming scheme || ) . 

2 The abbreviation in the title is derived from its authors' names. 



12 



We begin with an appropriate modification (rescaling) of the polynomial f{x). Let 

z,w = *?&<*>• (20) 

where x* = £*n is an integer parameter, < x* < n, and t = |_ nr o(£*)J • This choice is motivated 
by the following argument. The polynomial has to satisfy the inequality (14); any reasonable 
choice of Z(x) suggests that it be equal to h(x) at least at one point x = x*. This point is left 
a free parameter, chosen later. 



The program (20) gives rise to the following bound. For the reasons revealed in Q and 
outlined in footnote 3 below, in the quantum case the bound is valid for all but very low rates. 
Below ot\{q) is a certain small positive number dependent only on q. It can be computed for 
any q; for instance, ai(4) ~ 0.0028. We could not find a closed-form expression for it. 

Theorem 12 Let a x {q) < Rq < 1, R = (1/2)(1 - Rq). Then 



E(Rq, P ) 



< 



R - H q (5 l q p (l - R)) + T q (6 l q p (l - R),p) ai (q) < Rq < 2R l q p (p) - 1; 

R 2R l q p (p) -1<R Q <1. (21) 



Proof. We first prove the bound (|2l| ) and then prove feasibility of the program (P0[). By (45) 
we obtain 



j 2t h(x*) /n\ 2 /t + l + 7(n- t)\2 
aK?(q;x*)\t) V t+1 



z t (o ) = z y^^c) ( — ;,r » ) (22) 



Further, by 



q^hjx*) fn\ 

We would like to substitute these values in (|l6|). Recall the notation R 1 - = (1/n) log„ \C \ = 
(1/2)(1 + R Q ). Note that whenever 



log, {^)-f\cH 

i.e., H q (T (C)) +R X > 2fl,(Tb(r)) or 

R ± > H q (r (e)) 



we have Zt(0) = of Zq\C | I. The restriction R > iJ ? (ro(C*)) by our choice of r is equivalent 

to R 1 - > R q (£,*), which is always the case whenever £* > Sg P (R ± ). In this case the main term 
of the estimate (|lq) is given by the exponent of Zq\C |. Differentiating (l/ra)log„2:o on ^*, we 
obtain 

gg (i-p)e 
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The zero of this expression is £* = p, and d(\og zo)/d£* it is negative for £* G (p, 1). Thus, the 
optimal choice of £* is £* = p if this value is not less than 5 q {R) and £* = 5 q (R) otherwise. 



By (23), fl44|), and the fact that the exponent of h(x) is given by — T g (£,p), we obtain 

X - log, z = H,(r) - T q (C,p) ~ 1 - H q {r) + F g (f ) (24) 

Substituting this in ([jl]), we obtain (pl|). 

Let us prove that polynomial ( P0|) is admissible with respect to the restrictions (14)-(15). 
The proof will be broken into two cases, 

(a) £G (0,T] and 

(bHe(r,l). 

We begin with the first case and (14). We are only going to prove that it holds asymptotically, 
i.e., to prove the inequality 

^\og q Z t {in)<Uog q h(in). (25) 

Here we employ a method used in the corresponding part of Q. Namely, by our choice of r, 
the smallest zero (^2|) of Kt tends to £*; hence in the interval considered the exponent of Kt is 
given by (43). Then we can write 

-log q Z(Cn) = -T q (C,p) 

n 



/"£ (1 - y)j + y - qr + ^/((l - y)j + y - qr) 2 - 4jy(l - y) 

4 ° gq 2 7 (l-y) V - 

with o(l) terms omitted. Let ip(r, £) := 2 J*' * . . . dy + T q (^,p); then we have 

- log, Z(£n) + r,(e,p) = V(r,0 " </>(r,0- 
n 

Since (1/n) log„ Zt(x*) = (1/n) log„ h(x*), all we need to prove is that 



V (r, o =2 log (i-07 + e-gr + /(Mh+FrFW 3 !) 



27(1-0 

p_ 
(i-ph 



ipgg /n ^.. >o. o<e<r- 



First note that ip'Ar, £) is a monotone decreasing function of r = To(£*) and To (z) is a monotone 
decreasing function of z (|l9|). Hence if we prove that V4 is positive for £* = p this will also 

imply that it is positive for £* = ^(i? -1 ) > p. Therefore, let £* = p. In y, Appendix B] a similar 
function was proved to be positive. The proof proceeds as follows: consider the difference 



f c , [(i - £h + £ - qr + V((! - 0t + i - qr) 2 - 4 7 £(i - 



2 

p 



47 2 (l-0 2 (1~P)7' 
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£ 6 [0,£*),r = ro(p). The required inequality ipUroip),^) > is implied by g(To(p),S t ) > 0; 
the latter follows by the fact that g(To(p),p) = and that upon substituting r and simplifying 
we obtain a fraction whose denominator is positive and the derivative of the numerator on £ is 
negative in the whole segment £ £ [0, £*]. 

Now let us prove (15) in case (a). According to ||| the function ^[log Z t (£,n) — logz^ n ], £ £ 
(0, £*), achieves its maximum at £ = for 1 > Rq > oti(q)^. If p < 5lpi{R ± ) then as said 
above, we put £* = 5lpi(R ± ). This means that 

log q ^^H q ( t -)=R^e) = R\ 

and so -£o|C | — Z t (0) — » 0. Therefore, for any integer s, < s < x*, we have z s |C | — Z t (s) > 0. 

Finally if p = £* > 5 q p (R ) then (1/n) log g (Z((0)/zo) converges to a number less than i? . 
Hence for sufficiently large n and any integer s G (0, x*] we have ^IC -1 ! — Z t (s) > 0. This takes 
care of case (a) . 

To verify feasibility in case (b), i.e., to prove (14)-(15) for x* < x < n, we recall that a is 
the smallest zero of K t (q; x) + K t +i(q; x). Let y s be the smallest zero of K s (q; x). Then by the 



well-known properties of Krawtchouk polynomials we obtain that yt+i < a < yt\ so by (|42|), 
a \ y t +i as n — > cxd. Then we have, for large n and all integer x, x* < x < n, 



" l-\ „„\n—x. 



MX) <0<^;(1 -pT 

hence (14). To prove (15), observe that if x > x* then Zt(x) < and so for any integer s > x 
we have ZsIC" 1 "! — Z t {s) > 0. This exhausts case (b) and completes the proof. □ 

4.2 Hamming-type bound 

In the R-5 problem for nonbinary codes the bound [lj is not the best one known. It can be 
improved in several ways, in particular, in the frame of the polynomial method a better result is 
given in || . However, the technique in |2[ does not readily carry over to the present situation. 
Another, somewhat simpler bound that improves upon [|[ is the Hamming one which is better 
for a certain segment of rates close to 1. Therefore, in this subsection we derive a Hamming- 
type bound on E(Rq,p). This bound is valid for low error probabilities: p £ [0,p cr ], where 
the critical value p cr depends on q (it is 0.19 for q = 2 and 0.30 for p = 4). This improves 
Theorem [l^ for some values of Rq dependent on q and further extends the segment in which 
the exponent E(Rq,p) is known exactly. 



We begin with the polynomial 

n 

F e{x) =^2fiKi(q;x), 



i=0 



3 This is the reason for the bound to fail for very low rates both in [0 and here: the maximum shifts away 
from and the analysis becomes unmanageable. 
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where ft = K e (q;i) 2 . This polynomial is used in the proof of the Hamming bound on the size 
of the code with a given minimum distance ||. Our first goal is to show how to modify it for 
use in our problem. 

Delsarte |], p. 13] proved that 

n 
K i (q;x)K j (q;x) = ^2,p^K k (q; x) . 

k=0 

where p\, are the intersection numbers of the q-ary Hamming scheme. By a straightforward 
generalization of the binary case fi~2| , (A. 19)] we have 

v% = e L , , k . .) (" : ") ( 2k t 2a : ': j ) <« - »>«-*-*(« - d' 



s=0 



2k + 2s - i - j 



k + s — j 



Therefore, 



h = Y J PuKkW,i). 



k=0 

Substituting in F e (x), we obtain the following: 

n n n—k 



j=0 fe=0 s=0 



k \(n-k\ (2k + 2s-2e 

2k + 2s-2e)\ s )\ k + s-e 

.(q-2) 2e - 2s - k (q-iyK k (q;j)K 3 ( q -x) 

n n—k 

= EE 

fe=0 s=0 
■(q-2) 



k \(n-k\ (2k + 2s - 2e 

2k + 2s-2e)\ s )\ k + s-e 

n 

2e-2s-k 



\q-lYY,K k {q;j)K 3 {q; 
j=0 



X 



e-x/2 

E 



2x + 2s - 2e 



n — x\ (2x + 2s — 2e 

x + s — e 



(q-2 



,2e-2s-x 



(Q ~ I) 4 



s=max{0,e-i} 

where in the last step we made use of (40). 

Let us analyze the asymptotics of F e (x). Letting x = £n, s = <rn, and e = rn, we can write 
the exponent of the summation term as follows: 



n 



lo g 



n — x\ (2x + 2s — 2e 



2x + 2s — 2eJ \ s J \ x + s — e 
2£ + 2cr - 2r 



(?-2) 



2e-2s-x 



(9-1) 



£H 2 



i 



+ 2£ + 2o-2t 



\og q 2 + {l-i)H q 



l-£ 



+ ( 2r _2cr-£)log g (< ? -2)+C^~J. (27) 
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Computing the derivative of the last expression on a and equating it to 0, we arrive at the 
following condition: 

(2a + i-2rf(l-i-a){ q -l) _ 1 

a(2a + 2^-2r) 2 (q-2) 2 4' l ; 

It is not difficult to check (see pi, Appendix] ) that this equation has only one root in the interval 
max{0,r — £} < a < r — £/2. Denote this root by a®. The main term in F e {x) asymptotically 
corresponds to the value s = \jjqti\ • Thus, defining 



<t>(r,a,0 := I + (1 - 0H q [^-) + [^^^ =_j + 2 £ + 2o - 2rj log g 2 

+ (2r-2a-e)log g (g-2), (29) 

we observe that 

-log F e (x) = <j>(T,a ,£) + o(-). 

The analysis is complicated by the fact that gq itself is a function of r and £. 

Our general plan is, starting with F e (x), to construct a polynomial Z(x) so that Z(x) be 
equal to h(x) at one point and less than h(x) at all other integer points of the interval, thus 
guaranteeing feasibility. Together with (j^) this gives two conditions on the 2 parameters, 
o"o and r, both functions of £. It remains to make a suitable choice for £; this we simply guess 
prompted by an analogy in the binary case. This is the actual sequence of steps that we perform 
to derive the bound. Calculations, though elementary, are fairly involved, and we will not write 
them out in full. Instead, we perform a similar analysis in the binary case; this can be done 
explicitly within reasonable space and fixes ideas for the general result. 

THE binary CASE. We have the following simple expression for F e (x): 

x \ ( n — x 



(30) 



Fe{Xj - \x/2)\e-x/2 



£*ME(;_-^( t ^o- 



j=0 fc=0 

Note that F e {x) = when x > 2e. The exponents of F e (x) and h(x) are 

0(r,O:=ilog 2 F e (x) = l + e + (l-O^2(^|)+o(l), (31) 

-\og 2 h(x) = -T 2 (ti,p). 
n 

Now let x* = £*n be a point at which these exponents have equal slopes. Let us first convince 
ourselves that such a point exists and is unique. Indeed, the polynomial (J>'At,£) + (T2(£,p))'e 
is quadratic in £; its zeros are 

(l-p)(l-2r) 

VT=2p- 
Of them the one with the + sign is greater than one; the other one is always between and 1. 
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The equation [<^(t,£) = — (?2(^,p))t] f=f , defines r as a function of £*. Namely, 



r = r(p) = i(l-4^(l-n)- (32) 

Now we can define the polynomial Z{x) by rescaling F e {x) as follows. Let 

where e = [ttt,J . Note that we have ensured that Z{x) equals h(x) at x* and that their exponents 
are tangent; it will be seen below that for large n, Z(x) < h(x) at all the other integer points 
of the interval. It remains to choose the value of £*. This point is taken to maximize Zq, i.e., 



the estimate (16). Namely, the logarithm of zq equals 



- log 2 z = 2H 2 {t) - T 2 (e, P ) - <Kt, r ) + o(l). (33) 

n 

Substituting r from (^) and taking the derivative on £*, we find that the optimal choice is 
£* = p. Note that substituting £* = p and r from (^) in p3]), we find (1/n) log 2 zq — ► — 1. 

Let us examine feasibility of Z{x). As a preliminary remark, note that we are allowed to 
put £* = p as long as p is greater than the Hamming distance and we have to take £* equal to 
this distance otherwise. Indeed, substituting x = in (30), we obtain 

-log 2 F Tn (0)~l + fl- 2 (r) 

n 

\ log 2 h = \ ^g ^L-mJ ( 2 ; °) ~ 2F 2(r), 
the latter by fl45|) . From this and the definition of Z(x) it follows that whenever 

H 2 {t) > 1 - R ± , (34) 



we have Z(0) = o(zq\C \), needed for the estimate ( |lq ) to be nontrivial. So we can choose 
£* = p, r = r(p) if (p4j) holds and we choose r = il^ (1 ~~ -R -1- ) (the Hamming distance for the 
rate R- 1 ) and £* the root of i^CKC)) = 1 — -R 1 " otherwise. 



Note that by (|32| ) r grows on £*. Thus, (14) will follow in both cases if we prove that it 
holds for all (1/2)(1 — y/1 — 2p) < r < (1/2). As above, we are only going to prove that (14) 
holds asymptotically, i.e., that 

log 2 Z(x) - log 2 h(x) < 0, < $ < 2r, ^(1 - ^1 " 2p) < r < ^ (35) 

We begin with the case £* = p. Note that (1/n) log 2 /i(x) is a straight line; its derivative is 
log 2 (p/(l — p)) < 0. Inequality ( p5| ) will follow from the following set of conditions: 

(i) </£(t, £) = log 2 (p/(l — p)) has a unique zero for £ G [0, 2e]; 

(ii) o>^(r,0le=o>iog 2 (p/(i-p)); 

(iii) ^(r,OI ?=2e < log 2 (p/(l-p)) < 0, 
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where r = (1/2) (1 - yjl-2p). 

Condition (i) was established above. Conditions (i)-(iii) imply that 0|(r, £) < — T 2 (£,p) for 
£ 6 [0,p) and 0|(r, £) < — T 2 (£,p) for £ G [0,p). Indeed, suppose that 0(r, £) > log 2 (p/(l — p) in 
the neighborhood of £*, say on the left of it. This implies that in the small neighborhood of the 
point of tangency the derivative (pi is smaller that log 2 (p/(l — p)); however by (ii) it is greater 
than that for £ = 0, hence there is another point between and p at which (p and —T 2 have 
equal slopes, but this violates (i). Supposing that (p and —T 2 intersect at some point between 
and £*, we again find a similar contradiction. The second part of the claim follows by the same 



argument. Thus to establish (|35D it suffices to prove (ii)-(iii). 

We have 

4(t,£) = -log 2 ^ S^_ «; 

so 

4((1/2)(1 - - 2p, 0) = ^ log 2 2p > log 2 ^- (0 < p < i). 

This proves (ii). Condition (iii) is equally elementary; we omit the easy check. This establishes 
(14) for £* = p. 

Now suppose that £* > p. Condition (i) was proved above for any r. To prove (ii) and 
(iii) we only have to show that <pi (t, 0) grows and (pi (r, It) falls on r. Observe that (pi (r, 0) = 
\ log 2 4r(l — r) indeed grows as long as r < 1/2, which is true, and (pi(r, £) falls indefinitely as 
£ — > 2r. This proves (14), or rather (|35|), for q = 2. 



To prove (15), we choose the following tactics. We already know by (|34j ) that (15) holds for 
£ — ► 0. Hence it suffices to prove that the expression log 2 Z [\^n\) / z\^ n \ achieves its maximum 
at £ = 0. Numerical computations show that this condition holds if < r < t\ ~ 0.1069. 

Finally, by the definition, Z(x) = for x > [2rn]. Hence for these values of x (15) holds 
trivially. Otherwise by the preceding paragraph, (15) is true at least as long as i is less than 
the smallest root of K^ Tn ^ since otherwise the coefficients z\ = (K^ Tn ^(i)) 2 of Z(x) can be very 
small. The smallest zero is given by (|l9|); so the discussed constraint is satisfied in particular if 
r < t 2 , where r 2 = 0.1 is a root of 2r = (1/2) — yV(l — r). Thus a sufficient condition for (15) 
to hold true is that r < r cr = min{ri,r 2 } = 0.1. Note also that r{p) is monotone increasing in 
p. Hence we substitute £* = p in (|32| ) and denote by p cr the root of r(p) = r cr , p cr = 0.18. 



In summary, we obtain the following theorem. 

Theorem 13 Let R = (1/2)(1 - Rq). Then for any < p < p cr 

E(R )< [-l-R + T 2 {i,p) + ( p{H-\R),i) 1 - 2H 2 (r cr ) < R Q < 1 - 2H 2 (t( P )); 
Q,P) -{R 1-2H 2 (t( P ))<R q <1, 

where £ is a root of H 2 {t{^)) = R, and <p{-) is given by (31). 
This completes the argument in the binary case. 
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the general CASE. The analysis is similar but significantly more complicated since apart 
from r and £ we have a third parameter, a. In this part we are more sketchy than above. The 
polynomial Z(x) is again taken in the form 

„. , h(x*) „ , , 

where this time h(x) is as in (|lq) and x* = £*n is a parameter. 
By (29) 

(1/n) log g Z(x) = 0(t, ct , - T,(r ,p) - 0(t, ct , £*)• 

We proceed exactly as above. Namely, from the equation (f>Ur, 0"o>£*) = 0^g(£*>P))f we find t 
as a function of £* and o"o. This gives 

r = r «'^ := p(e-D(g-2) 



f ( g - i)(i - p)e 2 [(i - e) 2 (? - 1 -p) - o-o(i - e)(2g - 2 - «p) + ^(g - i)a - p) 

+ ~ p(e-i)(g-2) 

Next, we substitute this value of r in (|28| ) and find ctq as a function of £*. This gives 



^ 2g-2-gp-2(l-r)\/(g-l)(g-l-gp) ,„. 

^o := <r (£ ) = <? 2 (l-p) " ( } 

To complete the definition of the parameters we have to chose £*. As above, we take £* = p 
as long as this does not violate the feasibility condition (15)0. Substituting x = in (26), we 
obtain 

-log„F T „(0)~l + i? g (r) 

- log, /o = - log ^f rnJ (</; 0) ~ 2H q (r), 

n H n L J 

the latter by ((45|). From this and the definition of Z(x) it follows that whenever 

H q {r) >1-R L (37) 

we have Z(0) = o(zo\C ± \). 

So the best possible choice is £* = p, r = r(p, ctq) if (J37|) holds and r = H~ x {\ — R L ) 
(the Hamming distance for the rate R ) and £* the root of iJ ? (r(£, o"o)) = 1 — -R otherwise. 
Computations with Maple show that in the first case ^ log, zo — ► — 1, exactly as in the binary 
case above. We did not find a closed-form expression for the second case. 

Similarly to the binary case r should satisfy the inequality 

2r < 1 - Tll T _ -yVa-r) (cf. ©), 
g g g 



Though we do not prove this, this choice of £* is optimal with respect to the bound (|16| 
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or 



t < T2 :~- 



(3^q-2V2)(q-l) 



^q(9q - 8) 
Also similarly to the binary case we have to choose r such that 

r, Z(\£n\), , s 

arg max { log 2 — JJ } = 0. (38) 

For r = this maximum is obviously achieved at £ = (note that Fq(x) = q n 5 X Q and /^j = 1). 
Define ri,0 < n < T2, as follows: 

n = inf(r : ( |38[ ) does not hold). 

Let T cr = T\ if ri is well-defined and r cr = t 2 otherwise. The function again r{p) is increasing 
in p. Let p cr be the root of r(p, <Jo{p)) = T cr . Now we are ready to formulate the theorem. 



Theorem 14 Let R 



[1/2) (1 — Rq). Then for any < p < p cr 

-1-R+ T q (C, P ) + (l>(H- 1 (R),a (e),^ 



E(R Q ,P) < { 



1 - 2H q (T cr ) <R Q <1- 2H q (T(p,o- )); 
R, l-2H q (T(p,a ))<RQ<l, 

R, o"o is defined in (|36|), and 4>{-) is given by (29). 



where £* is a root of Hq(r(£, 0o(£))) 

Remark For q = 4, p cr = 0.301, and numerical computations show that ( p8| ) holds true in 
the entire interval r S [0, t 2 ]. Therefore in this case r cr = r 2 . 

Fore reference purposes we composed a short table of values of the bounds for q = 4, p = 0.1. 



Rq 


Existence 


A-MRRW 


Hamming 





0.5260 


0.6270 


- 


0.1 


0.4637 


0.5458 


- 


0.2 


0.4054 


0.4685 


0.4774 


0.3 


0.3509 


0.3952 


0.3951 


0.4 


0.3 


0.3262 


0.3216 


0.5 


0.25 


0.2618 


0.2567 


0.6 


0.2 


0.2028 


0.2003 


0.7 


0.15 


0.15 


0.15 


0.8 


0.1 


0.1 


0.1 


0.9 


0.05 


0.05 


0.05 


1 












These bounds are also plotted in Fig. 1. It can be seen that the Hamming bound is the 
best of the two upper bounds for large rates. Unlike the classical case, the upper bounds do 
not approach the lower bound as the rate Rq becomes small. However this is due rather to the 
way of measuring the rate of quantum codes than to an imperfection of the method. Indeed, 
roughly speaking, the case Rq = corresponds to classical codes of rate R = 1/2 (cf. (§)). The 



function E(Rq,p) is known exactly at least for 2R£(p) — 1 < Rq < 1. In fact, by Theorem 14 



lp. 



the left end of this interval is provably smaller than this value; however, it is difficult to make 
any exact statements other than just plotting the bounds. 
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5 Appendix 

Let Ki(q;x) be the g-ary Krawtchouk polynomial, 7 = q — 1. Here we list its properties used 
in the paper. 

The following 3 basic facts are well known: 

n 

Y^ Ki(q; x)z l = (1 + jz) n ~ x (l - zf (x integer); (39) 

i=0 

n 

Y J K r {q;i)K i {q-s)=q n 5 rs - (40) 

i=0 
t n 

f{x) = Y J UKi{q]x) «> /^g-^/Ol^fei), (41) 

i=0 j=0 

where in (41) f(x) is any polynomial with deg/ < n. 

Let y s be the smallest zero of K s . For s = an, n — > 00 we have |j]], 

-=T(a)+o(l), (42) 

n 

where the function r(-) is defined in (|J9|). Further, by @ for n — > 00 and £ E [0,t(<t)] 

1 



log# s (g;£n) ~ HJa) 



n 

, /" g w C 1 - z/h + 2/ - gg" + V((i - s/b + y - g^) 2 - 4 7y(! - r y) d ( 43 \ 

Jo 2 7 (1 -y) 

In particular, for c = r(£), i.e., £ = r(<r), 

1 -lo gq K, n ( q ^n) = 1 + H ^ ) 2 - H ^ ) + o(l). (44) 

Finally, from the definition of K s we find 

K s (q;0)=(^\ 7 s . (45) 
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Figure 1: Bounds on E(Rq,p) : Existence bound, Theorem [j|; Aaltonen-MRRW-type bound, 



Theorem 12; Hamming- type bound, Theorem 14; p = 0.1. 
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